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ABSTRACT 

Using the archive data of the H 13 CO + (J = 1 — 0) line emission taken with the 
Nobeyama 45 m radio telescope with a spatial resolution of ~ O.Olpc, we have identified 
68 dense cores in the central dense region of the p Ophiuchi main cloud. The H 13 CO + 
data also indicates that the fractional abundance of H 13 CO + relative to H2 is roughly 
inversely proportional to the square root of the H2 column density with a mean of 
1.72 x 10 -11 . The mean radius, FWHM line width, and LTE mass of the identified cores 
are estimated to be 0.045 ± 0.011 pc, 0.49 ± 0.14 km s -1 , and 3.4 ± 3.6 M , respectively. 
The majority of the identified cores have subsonic internal motions. The virial ratio, 
the ratio of the virial mass to the LTE mass, tends to decrease with increasing the LTE 
mass and about 60 percent of the cores have virial ratios smaller than 2, indicating that 
these cores are not transient structures but self-gravitating. The detailed virial analysis 
suggests that the surface pressure often dominates over the self-gravity and thus plays a 
crucial role in regulating core formation and evolution. By comparing the p Oph cores 
with those in the Orion A molecular cloud observed with the same telescope, we found 
that the statistical properties of the core physical quantities are similar between the two 
clouds if the effect of the different spatial resolutions is corrected. The line widths of 
the p Oph cores appear to be nearly independent of the core radii over the range of 0.01 
— 0.1 pc and deviate upwards from the Heyer & Brunt relation. This may be evidence 
that turbulent motions are driven by protostellar outflows in the cluster environment. 

Subject headings: ISM: clouds — ISM: individual (p Ophiuchi) — ISM: structure - 
stars: formation — turbulence 



1. Introduction 

Observation s of young stellar populations in the G alaxy have reve aled that th e majo rity of stars 



form in clusters (jLada &: Ladall2003l ; lAUen et al.ll2006l ). For example. lLada et al.l (|199ll ) performed 
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an extensive near-infrared imaging survey of the central regions of LI 630 in Orion B and found 
that about 60 — 90 % of all the young stellar populations associated with the mo lecular cloud form 
only in three rich embedded clusters. Using the 2MASS point source catalog, ICarpenterl (j2000) 
estimated the fraction of young stellar populations contained in clusters to be 50 — 100 % for 
nearby cluster forming molecular clouds such as Perseus, Orion A, Orion B, and MonR2. Recent 
near-infrared surveys of young stellar populations using the Spitzer Space Telescope hav e confirmed 



Poulton et al 



2008 



that clustered star formation is the dominant mode of star formation in the Galaxy (lAllen et al 
20061: 



Roman-Zuniga et al.ll2008l ). However, previous studies of star formation 



have focused on star formation in relative isolation, for which the environmental effects that play 
a crucial role in clustered star formation are thought to be minor or at least secondary. Thus, how 
stars form in the cluster environment remains only poorly understood. 

An important clue to understanding star formation in clusters has come from recent millimeter 
and submillimeter observations of the nearby cluster-forming region, the p Ophiuchi molecular 
cloud (L1688), wh i ch uncovered dense cores with a mass spec trum that resembles the Salpeter IMF 



(Motteetal 



1998; 



Johnstone et al 



2000a 



Stanke et al.ll2006l ) . Similar studies in other star forming 
regions have confirmed that such a resemblance between the core mass spectrum and the stellar 



IMF seems to be common in nearby star for ming regions (ITesti & Sargentlll998l : lMotte et al.l 



Reid k Wilson 



2006 



Ikeda et al. 



2007 



2001 



20091 ) . These observations sug gest that the stella r IMF may 



be determined to a large extent by the core mass distribution (e.g., I Andre et al.l 120071 ). although 
mass accretion from ambient gas onto forming stars may s ometimes plays an important role i n 
determining the final masses of stars (jBonnell et al.1 l200ll : iBate et al.l 120031 : IWang et al.l 120091 ) . 
Therefore, it is important to investigate the formation and evolution of dense cores in the cluster 
environment. 

Recent theoretical studie s on core formation in the cluster environment have emphasized the 



2005 



role of supersonic tur bulence (IKlessen et al.lll998l ; iTillev h Pudritz 



2004 



Vazquez-Semadeni et al 



Dib et al.l 120071 ) . These studies have demonstrated that dense cores are formed in regions 



compressed by converging turbulent flows. However, the origin of this supersonic turbulence and 
its role in core formation are still a matter of debate. Two main scenarios have been proposed. In 
the first scenario, formation of dense cores (and thus star formation) is considered to be completed 
only in one turbulent crossing time, i.e., o ne dynamical time, so that the cluster forming region s 
are destroyed within one dynamical time (lElmegreenl |2000j; iHartmannl l200ll : iHeitsch et al.l 120081 ) . 
In this case, the mass spectrum of de nse cores is shaped by the cascade of primordial interstellar 
turbulence that has a power law form (|Padoan fc NordlundlEooi ') . On the other hand, in the second 
scenario, star formation is considered to continue at least over several dynamical times. In this case, 
additional turbulent motions must be replenished to keep supersonic tu rbulence in molecular clouds 
becau se supersonic turbulence dissipates quickly in one dynamical time ([Stone et al.lll998l ; lMac Low 
19991 ). In such circumstances, dense cores are continuously created (and destroyed) in molecular gas 
disturbed strongly by the additional turbulent motions. In the cluster environment, stellar feedback 
such as stellar winds, HII regions, and protostellar outflows has been discussed as the dominant 



-3 - 



sourc e of additional supersonic turbulence (jNorman fc Silkl I198Q ; iMatznerl 120071 ; iLi &: Nakamura 
2006). Other sources of turbulence such as supernova bla stwaves, external cloud shearing, or 
converging flows of cloud formation ma y also play a rol e (e.g., iMcKee fc Ostrikerl 120071 ). As for the 
protostellar outflow-driven turbulence, ISwift fc Welch! (|2008l ) identified the characteristic energy 
injection scale at ~ 0.05 pc for the L1551 star forming region where a group of stars are forming, 
indicating that t he protostellar outf lows are the primary source of regenerating turbulence in this 
region. Recently, lOffner et al.l (|2008l ) performed three-dimensional (3D) hydrodynamic simulations 
of turbulent molecular clouds and suggested that the physical properties of dense cores may be 
affected directly by the presence or absence of turbulent feedback, although further justification 
may be needed. Therefore, to shed light on the issue of core formation in the cluster environment, it 
is essential to investigate the kinematics of dense cores. In this paper, we thus analyze the physical 
properties of dense cores in the nearest cluster forming region, the p Ophiuchi molecular cloud. 

The p Qphiuchi molecula r cloud is the nearest cluster forming regi ons at a distance of ~ 125 



pc (e.g., lLombardi et al 



2008 



Loinard et al. 



2008 



Wilking et al J 12003 ). The dense part of the p 



Ophiuchi molecular cloud (hereafter the p Ophiuchi main cloud) is known to harbor a rich cluster 



of young stellar objects (Y SOs) in different evolutionary stages (jWilking et al J 120081 : Enoch et al 



2008 



J0rgensen et al.ll2008l ). Recent millimeter and submillimeter dust continuum observations of 



the p Ophiuchi main cloud have revealed that a large number of dense cores are concentrated in 
the central dense part of the cloud and th eir mass spectra are similar in shape to the stellar IMF 



(jMotte et alJll998l ; |johnstone et al 



2000 aJ) . By comparing w ith isothermal, pressure-confined, self- 



Johnstone et al.1 (|2000al ) argued that the dense cores detected 



gravitating Bonnor-Ebert spheres, 
by t he dust continuum em ission are likely to be self-gravitating and expected to form stars (see 
also I Johnstone et al.1 12004 ) . However, to accurately assess the dynamical states of the dense cores, 
it is necessary to perform molecular line observations using high density tracers such as H 13 CO + 
( J = 1 — 0) and N2H + ( J = 1 — 0), both of which have the critical densities for excitation as high 
as 10 5 cm -3 . Recently, the Nobeyama observatory has released the archive data of the H 13 CO + 
( J = 1 — 0) line emission toward the p Ophiuchi main cloud taken with the 45 m radio telescope. 
At a distance of 125 pc, the spatial resolution of the data of 21", corresponding to ~ 0.01 pc, is 
comparable to those of the dust continuum observations, e.g., the SCUBA at the JCMT with a 
beam size 14" at 850 pm. Thus, the Nobeyama archive data is suitable for analyzing the dynamics 
of small scale structures comparable t o the ones det e cted by the dust continuum observations. Our 
analysis is complementary to that of lAndre et al.l (|2007l ) who observed several subclumps (Oph 
A, Bl, B2, C, E, and F) with the IRAM 30 m radio telescope us i ng No H+ (J = 1 - 0) line to 
measure the line widths of the dust cores detected bv Motte et all (H998I ). The H 13 CO+ data of 



the Nobeyama 45 m telescope was taken in a larger area cove ring the who l e den se region of the p 
Ophiuchi main cloud, including all the subclumps observed bv I Andre et al.l (|2007l ). Since the dense 
gas detected by molecular line emissions such as H 13 CO + ( J = 1 — 0) and N2H + ( J = 1 — 0) may 
not be alway s well associated w ith the dust cores and therefore the line widths of the dust cores 
estimated bv lAndre et al.l (|2007l ) may not be accurate, we identified the dense cores directly from 
the molecular line data cube and estimated their physical quantities such as the core radius, mass, 
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and line width. 

The rest of the paper is organized as follows. First, the description of the archive data is 
presented in £}2j In ^3] we show the spatial distribution of the dense gas detected by the H 13 CO + 
emission, which indicates very clumpy structures. Then, in £JH we identify the clumpy structure s 



as dense cores applying a clump-finding algorithm, the so-called clumpfind ([Williams et al.lll994l ). 
and derive the physical properties of the identified cores. To accurately evaluate the core masses, 
the fractional abundances of H 13 CO + relative to H2 are derived directly by comparing the H 13 CO + 
column densities with the H2 column densities obtained from the 850 pm dust continuum emission. 
In £}5] we discuss the mass spectrum of the identified cores and their dynamical states. We also 
compare the physical properties of the identified cores with those in Orion A, taking into account 
the effect of different distances (or different spatial resolutions). Finally, we summarize our main 
conclusion in 36l 



2. Data 

We used the archive data of the p Ophiuchi main cloud taken in H 13 CO + (J = 1 — 0) molecular 
line (86.75433 GHz) with the Nobeyama 45 m radio telescope. The fits data is available from the 
web page of the Nobeyama Radio Observatory at the National Astronomical Observatory of Japan 
(http://www.nro.nao.ac.jp/). The observations were carried out in the period from 2002 March 
to 2003 May. All spectra were obtained in the position-switching mode. The data has 118 x 98 x 60 
grid points in the a-5-v^sR space, covering the whole dense region of the p Ophiuchi main cloud 
with size of about 1.5 pc xl.3 pc. The grid size of the data was 21", corresponding to 0.013 pc 
at a distance of 125 pc. The main beam efficiency of the telescope was 0.51. We note that the 
FWHM beam size of the telescope was 18" at 87 GHz, and therefore the data were almost full- 
beam sampled, considering the typical pointing error of a few arcsecond (see the web page of the 
Nobeyama Radio Observatory). The velocity resolution is 0.13 km s -1 and the average Root Mean 
Square (rms) noise determined from signal-free channels is a = 0.11 K in 



3. Distribution of Dense Molecular Gas in the p Ophiuchi molecular cloud 

3.1. Overall Distribution of the H 13 CO + emission 

Here, we present the overall distribution of dense gas observed by the H 13 CO + ( J = 1 — 0) 
emission. Figure [Ija, shows a total integrated intensity map of the H 13 CO + emission toward the p 
Ophiuchi main cloud. For comparison, the 850 pm. image obtained with the SCUBA at the JCMT 
is shown in Figure [T]o. While the distribution of the H 13 CO + emission appears to cover a larger 
area than that of the 850 ^m emission, the overall spatial distribution of the H 13 CO + emission 
is similar to that of the 850 pm emission. We note that it is very difficult to observe extended 
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structures with bolometers mounted on ground-based telescope due to the atmospheric emission. 
The extended structures larger than ~ 1.5' (~ 0.06 pc) are therefore suppressed in removing the 
atmospheric emissi on during the data reduc tion process, making the 850 /um map mostly devoid of 
extended emission (j Johnstone et al.ll2000bl ). Our H 13 CO + map shows that the well-known dense 
subclumps (Oph A, Bl, B2, C, E, and F. Note that the Oph D region is out of our area.), previously 
identified by D CO + and other ob s ervations, are very clum py as recent submillimeter observations 
have revealed (jMotte et al.l Il998l : I Johnstone et al.1 l2000al ). In the next section, we identify the 



clum py structures as dense cores applying a clump-finding algorithm, clumpfind (jWilliams et al 



1994 ). and derive physical properties of the identified cores. 



3.2. Individual Regions 



Blow-up total integrated intensity maps toward the Oph A and Oph Bl, B2 and C regions 
are shown in Figures (2£t and [2b, respectively. Iso antenna temperature surfaces of the same areas 
in the 3D data sp ace («-<5-ut,rr ) are also shown in FigureEl Our intensity maps can be compared 
with Figure 2 of lAndre et al.1 (|2007l ). the N2H + (J = 1 — 0) intensity maps obtained with the 
IRAM 30 m telescope. The beam sizes for both the telescopes are comparable: the HPBW of 
the IRAM 30 m, 26.4", is about 1.5 times as large as that of the Nobeyama 45 m telescope, 18". 
Both the molecules are also expected to trace dense molecular gas in nearly the same density range 
because the critical densities for excitation of H 13 CO + ( J = 1 — 0) and N2H" 1 " ( J = 1 — 0) lines are 
comparable (?^cr,H 13 co+ = 8 x 10 cm -3 and re cr n 2 h+ = 2 x 10 5 cm -3 ). In fact, the distribution of 
the H 13 CO + emission is remarkably similar to that of the N2H + emission. 

The most intense line emission of H 13 CO + (J = 1 — 0) comes from the Oph A region. The 
maximum antenna temperature r eaches T*^ ~ 4 K a t the position that coincides with that of a 
starless object, SMI, identified by I Andre et al.l (|l993l ) within our spatial resolution. A prominent 
feature of this region is a dense ridge running from north to sout h, which is a l so see n in the 850 
//m map. The starless objects SMI, SM1N, and SM2 identified bv lAndre et al.l (|1993l ) are located 
just at the densest part of the ridge. Although relatively strong H 13 CO + emission is detected 
toward VLA1623, the prototypical Class object, it is difficult to distinguish the H 13 CO + emission 
associated with VLA 1623 from the ambient component because of our resolution as large as 21" . 



In the 850 um map , IWilson et al.l (Il999l ) found two filaments in the Oph A region [see Figure 



1 of I Wilson et al.1 ( 19991 )] . Similar structures can be recognized in our H 13 CO + map. In the 850 
jj,m map, the two filaments (the dashed lines in Figure [2h.) intersect at the position (RA, Dec) = 
(16 h 26 m 29 s ,-24°22 / 45 // ; J2000). In the H 13 CO+ map, the north-east filament is not seen, whereas 
the north-west filament, a part of the Op h A ridge, has strong H 13 C O + emission. Such a feature is 
consistent with that of the N2H + map bv lDi Francesco et al.l (|2004l ). indicating that the north-east 
filament presumably consist of diffuse warm gas. The filament is reminiscent of the Orion bar where 
the dust continuum emission enhanced by heating from PDRs and young stars is prominent but 
the H 13 CO + line emission is not seen ( Ikeda et al.1 120071 ). In fact, the location of the north-east 
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filament agrees well with the edge of the PDR seen in ISO images ([Abergel et al.lll996l ) 



In the western side of the Oph A ridge. IWilson et al.l (| 19991 ) found another interesting structure: 
two arcs that are labeled with No. 1 and 2 in Figures [2k. and [3k.. In addition to these two, we 
identified 3 new arcs in the Oph A region (labeled with No. 3 through 5). Although the No. 1 arc 
is not clearly recognized as an arc in the total intensity map of Figure [2k. owing to the presence 
of the extended diffuse component, the No. 1 arc as well as the No. 2 and 3 can be clearly seen 
in the iso-temperature surfaces in Figure [3k- These arcs may be created by the stellar wind from 
the young B3 star SI. Another possibility is the effect of three outflows detected in this region 
whose o utflow axes are almost parallel in the a-5 space and inclination angles are around 60 — 80 
degrees (jKamazaki et al.ll2003l ). One of the three outflows is the most spectacular outflow from 
the Class object VLA 1623. Both sides of the No. 1 arc extended from VLA 1623, where no 
H 13 CO + emission is detected, appear to be parts of a cavity created by the VLA 1623 outflow. 
On the other hand, in the eastern side of the Oph A ridge, we can recognize two new arcs. These 
arcs (labeled with No. 4 and 5) appear to be related to a giant protostellar outflow that has been 
recently discovered by Nakamura et al. (2009, in preparation) on the basis of the 12 CO ( J = 3 — 2) 
observations. 

A couple of features that are likely to be created by protostellar outflows are also seen in 
the Oph B2 region, the densest region after the Oph A. The H 13 CO + map s hows a very clumpy 
structure with ~ 10 local peaks in the region. Many dust cores identified by iMotte et al.l (119981 ) 
are distributed mainly in the bright area of the H 13 CO + map, where some cores identified by 
the clumpfind in the data cube (see $3$ are overlapped on the plane of the sky. It is difficult to 
distinguish such an overlapped core from the other in the H 13 CO + total integrated intensity map 
or the dust continuum map that does not have velocity information. In the southern part of the 
Oph B2 region (labeled with No. 1), a hole is seen in both the H 13 CO + and 850 p,m maps. The 
hole is particularly clear in the data cube in Figure [3p. Two other holes are also seen near the 
northern (No. 2) and southern (No. 3) edges of the B2 region. These holes are likely to be created 
by the giant protostellar outflow recently found by Nakamura et al. (2009). The blue lobe of this 
outflow is extended beyond the B2 region, reaching the No. 4 and 5 arcs. The apparent luminosity 
and the total mass of this outflow are largest in p Oph. The driving source is likely to be a class 
I YSO Elias 32 or 33. No. 1 and 2 holes are likely to be created by the red and blue lobes, 
respectively. The interaction between the protostellar outflows and the dense gas will be discussed 
elsewhere (Nakamura et al. 2009, in preparation). Cavities that could be created by protostellar 
outflo ws are also found in NGC1333, a nearby cluster forming clump in the Perseus molecular 
cloud |Knee fe Sandelll boocl : ISandell & Kneel boool : buillen et aDbood ). suggesting that density 
structures in cluster forming regions are strongly affected by protostellar outflows. Nakamura et 
al. (2009) also found other examples of the dynamical interaction between the dense gas and the 
outflows in p Oph, which will be shown elsewhere. 
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4. Physical Properties of Dense Cores 



We identified 68 H 13 CO + cores in the p Ophiuchi main cloud, using the clumpfind method 
described in §4.11 Here, we describe physical properties of the identified cores. The definitions of 
the physical quantities of the cores are presented in §4.21 The fractional abundance of H 13 CO + 
relative to H2, which is important to the determination of core mass, is estimated in §4.31 In 
§4.4^ we describe the statistical properties of the physical quantities of the identified cores, and 
comp are their positions with those of YSOs identified by the Spitzer Space Telescope mid-infrared 
data (IJorgensen et al, l2008l ). We also compare our cores with 850 /im cores listed in Table 7 of 
Jorgensen et al.1 (120081 ) . Finally, we show in §4.61 the correlation among physical quantities such as 
the line width-radius and virial ratio-mass relations. 



4.1. Core Identification 



To identify dense cores f rom t he H 13 CO + data cube, we applied the clumpfind method de- 
scribed in §3.2 of llkeda et al.l (|2007l ). In brief, we adopte d the threshold ( TA,th ) and stepsize (AT^) 
of the 2 a noise level for core identification, following I Williams et al.l (119941 ) who recommended 
using this set of the parameters based on the results of their test simulations. To identify the cores 
with reasonable accuracy, we imposed the following 3 criteria: (1) a core must contain two or more 
continuous velocity channels, each of which has at least 3 pixels whose intensities are above the 
3a noise level, and (2) the pixels must be connected to one another in both the space and velocity 
domains. (3) Furthermore, we rejected core candidates located at the edge of the observed region. 
The total intensity of the minimum mass core identified by this method has about 60 a levels. Using 
our definition of core mass (eq.pQ), a minimum core mass is estimated to be 0.28 M & and 0.38 M Q 
for T ex = 12 K and 18 K, respectively. The total mass of the identified 68 cores is estimated to be 
228 Mq, 60 % of the total mass detected above the 2a level. 

In general, the physical properties of the identified cores depend on the core identification 
scheme and its parameters. In our case, we expect that our clumpfind can detect the real 3D 
density structures in the cloud from the H 13 CO + (J = 1 — 0) spe ctral line posit io n-pos ition- 
velocity data cube because of the following reasons. According to IWilliams et al.l (|1994l ). the 
physical properties of the cores identified by clumpfind are affected siginificantly by the adopted 
threshold and stepsize unless they are carefully chosen. They recommend using the 2 a noise level 
for both the threshold and stepsize, for which the clumpfind could detect and accurately measure 
the individual core properties in their test simulations. We adopt this set of the parameters in our 
clum p find. Recently , sever al authors revisted this problem (IPineda et al.1 120091 : llkeda &: Kitamura 
2009J). IPineda et al.1 (|2009l ) reexamined the effect of varying the values of the threshold (Ta^Ii =3, 5, 
and 7 a) and stepsize (3a < AT4 < 20a) for the core identification, and confirmed that the physical 
quantities of the identified cores depend on these parameters especially when the distribution of 
the emission is not localized in the data cube and the larger values of the parameters are used 
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[see also ISchneider fc Brooks! (|20Q4j ) for the case of the low-spatial resolut i on (~ 1 pc), optically 
thick 12 CO (J = 1 - 0) data with T A ,th = AT A = 3 a]. Ilkeda fe Kitamural feood ) explored more 
reasonable parameter ranges of T^th = 2 — 3a and 2a < AT A < 3a suitable for their core study 
in the OMC-1 region, and demonstrated that the power-law in dex of a core mass spectrum does 
not depend on the threshold and stepsize. I Williams et al.l (|1994i ) also discussed the "overlap" effect 
for which multiple cores that almost overlap along the line-of-sight are not separately identified 
by the clumpfind when they have nearly the same line-of-sight velocities. This overlap effect is 
likely to be significant when the gas is highly turbulent and the emission tends to be distribute d 
in the large v olume of the 3D data cube (jBallesteros-Paredes Mac Lowll2002l ; ISmith et al.ll2008l ). 
For example, lOstriker et al.l (|200ll ) performed 3D MHD turblent simulations and claimed that the 
emission lines are often singly peaked even in the presence of multiple condensations along the 
line of sight. We note that their model clouds are highly turbulent, gravitationally unbound and 
the condensations have extremely large virial ratios, ranging from 10 to 10 3 . This suggests that 
the overlap effect may be crucial for less dense gravitationally unbound molecular gas components 
observed by 12 CO and 13 CO (see lOssenkopf Mac Lowll2002l . for the case of the Polaris molecular 
cloud, a nearby diffuse, gravitationally unbound molecular cloud.). In contrast, for our H 13 CO + 
data, the overlap effect is expected to be minor because the H 13 CO + emission can trace only dense 
gas with densities greater than ~ 10 5 cm -3 and appears to be reasonably localized in the 3D data 
cube, indicating that the dense gas is likely to have a small volume filling factor in the cloud and 
moderately quies cent (the Ma ch number ~ 2). Note that the low density gas has the Mach number 
of about 5 — 6 (|Lorenl ll989). In fact, our cores have virial ratios of about 2, as shown in ^4.51 
Thus, we expect that the clumpfind can identify the 3D density structures reasonably well from 
our H 13 CO+ data. 



A caveat of our analysis is that the cores identified by the above procedure (and any other 
core identification schemes) are resolution dependent in the sense that the substructures smaller 
than the telescope beam size cannot be resolved. It remains unclear whether an identified core is 
a unit that is separated dynamically from the background media. However, we expect that the 
structures identified from high density tracers such as H 13 CO + (J = 1 — 0), N2H + (J = 1 — 0), 
and dust continuum data correspond to the real dense structures in the cloud as long as the spatial 
resolution is less than ~ 0.01 pc, because recent numerical simulations of pc-scale cluster forming 
clumps demonstrated that the dense structures like filaments often created i n the clump have sizes 



of order 0.1 pc or less along the minor axes and densities of order 10 5 cm 3 ( Li fc Nakamura 



2006 



Nakamura fc Lill2007l ). Therefore, in the following, we discuss the statistical properties of the cores 
identified by the above procedure, assuming that the identified cores can be regarded as units that 
are reasonably separated dynamically from the background media. Further investigation of this 
problem should be done when higher spatial resolution observations are available. 
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4.2. Derivation of Core Properties 



We derive physical properties of the cores using the definitions described in §3.3 of llkeda et al 



(|2007l ) . In brief, the position and local standard of rest (LSR) velocity of each core are determined 
at the pixel with the largest antenna temperature, T* A peak , within the core. The core radius, R CO re, is 
determined by taking the projected area enclosed by the 2 a level contour and computing the radius 
of the circle required to reproduce the area, taking into account the correction for the telescope 
beam (18"). The aspect ratio is computed as the ratio of the major to minor axis lengths that are 
determined by the two-dimensional Gaussian fitting to the total integrated intensity distribution of 
the core. The FWHM line width, dv COTe , is corrected for the velocity resolution of the spectrometers 
(= 0.13 km s _1 ) and for line broadening due to hyperfine splitting of H 13 CO + . The LTE mass is 
estimated as 



M LTE 



x 10 

D 
125pc 



-2 



X 



H 13 CCH 



1.72 x 10 

2 / A 0\ 2 

21" 



-ii 



-i 



0.51 



-i 



T ex exp (4.16/T t 



K km s 1 



M G 



(1) 



by assuming that the H 13 CO + emission is optically thin. Here, X-^wqq+ is the fractional abundance 
of H 13 CO + relative to H2, T ex is the excitation temperature, D is the distance to the p Ophiuchi 
main cloud, ^Tjj^Aui is the total integrated intensity of the core and Avi = 0.13 km s . We 
adopt the fractional abundance X H i3co+ = 

1.72 x IP" 11 (see &T3l). and t he distance D = 125 pc on 
the basis of the values r ecently updated by Lombardi et al. (feoosl ) and Loinard et alj ( 2008) (see 
also Wilking et al.ll2008l ). The excitation temperature is assumed to be T cx = 12 K (jMotte et al. 
19981 ) except for the area shown in Figure [2k, i.e., t he Oph A region, where T cx = 18 K following 
the N2H + observations by lDi Francesco et al.l (J2004J). 



The virial mass M v i r is calculated as 



5a 1 R ca redv 2 



tot 



209a 



In 2G 

Rc 



pc 



dv t , 



ot 



km s 1 



M , 



and 



dv 



tot 



dv r 



+ 8ln2k B T 



1 



1 



1/2 



(2) 



(3) 



K pm H m obs , 

where kg is the Boltzmann constant, T is the kinetic temperature of the molecular gas, fi is 
the mean molecular weight of 2.33, run is the mass of a hydrogen atom, m b s is the mass of a 
H 13 CO + molecule, and a is a dimensionless para meter of order unity whic h measures the effects 
of a nonuniform or nonspherical mass distribution (IBertoldi fc McKedll992l ). For a uniform sphere 
and a centrally-condensed sphere with p oc r -2 , a = 1 and 5/3, respectively. For our cores, the 
effects of the nonspherical mass distribution appear to be small because the aspect ratios are not 



- 10 - 



so far from unity (see also Figure 2 of Bertoldi & McKee 1992). The virial ratio, the ratio between 
the virial mass to the LTE mass, is defined as a v ir = -^vir/-^LTE) which is equal to twice the ratio 
of the internal kinetic energy to the gravitational energy. 

The mean number density of the core, n, is calculated as the LTE mass divided by the volume 
of a sphere with radius, R CO re, 

_ = 3M LTE 



4.3. Fractional Abundance of H 13 CO 



To evaluate the fractional abundance of H 13 CO + , we smoothed the 850 //m map with the 
angular resolution of 21" to match the H 13 CO + total integrated intensity map and computed the 
column densities of H 13 CO + (A r H i3co+) an d H2 (-/Vh 2 ) at each pixel in the maps. Here, we computed 
-/Vh 2 from the relation 

^H 2 = S%!P/[n hcamf im H K S50 B S50 (T)}, (5) 

where S^ m is the 850 [j,m flux density per beam, fibeam is the main beam solid angle, B$5o(T) 
is the Plan k function for a dust temperature T, and Kano (= 0.01 cm 2 g -1 ) is the d ust opacity 



at 850 /im (jHenning et al. 



1995 



Ossenkopf Henning 



1994 



Johnstone et al.1 l2000al ). The dust 



temperature is assumed to be equal to the gas temperature. In Figure H] the fractional abundances 
X H i3co+ computed at all the pixels above the 3 a levels for both the H 13 CO + and 850 p maps 
are plotted against Nn 2 • We note that the H2 column density derived from the 850 fim map is very 
sensitive to the temperature and the dust opacity assumed, and therefore have uncertainty of at 
least a factor of a few. The fractional abundance of H 13 CO + tends to decrease with increasing H2 
column density. The best-fit power-law is given by 

23, 



logX H i3 CO + = (-10.87 ±0.01) + (-0.514 ± 0.024) log(iVH 2 /10 cm _i! ) 



(6) 

— 1/2 

with a correlation coefficient 1Z = 0.68, indicating that X^{i-aqq+ is nearly proportional to iV H2 
This tendency is in good agreement with a theoretical consideration of the ionization fraction in 
molecular clouds if the column density is proportional to the local volume density for each molecule. 
When the ionization rate by cosmic rays (oc n n , where n n is the density of neutral gas) balances 
with the recombination rate (oc njn e , where rii and n e are the densities of positive ion and electron, 
respectively, and rii y n P ), the ionizat ion fraction is inversely proportional to the square root of the 
neutral gas density (|Elmegreenlll979l ). In typical regions in molecular clouds, the most abundant 
positive molecular ion is expected to be HCO + . Therefore, the fractional abundance of H 13 CO + as 
well as that of HCO + is expected to be inversely proportional to the square root of the neutral gas 
density if the fractional abundance of H 13 CO + relative to HCO + is almost constant. More detailed 

— 1/2 

studies indicate that the relation X h i3qo+ oc n n is a reasonable approximation of the fractional 
abundance in the density ran ge of n < 10 7 cm" 2 , although metal ions such as Mg + may be the 
most abundant positive ions (jNishi et al.lll99ll ; iNakano et al.ll2002l ). We note that as mentioned 
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in <j3l in the 850 pm data the structures larger than 1.5' (~ 0.06 pc) are removed during the data 
reduction process. To evaluate how the effect of this artificial filtering affects the estimate of the 
H 13 CO + abundance, we repeated the same analysis by filtering out all the structures larger than 
1.5' in the H 13 CO + map and confirmed that the effect of this artificial filtering is negligible. This 
is because the dense gas detected by the 850 pm and H 13 CO + is spatially well- localized. 



In nearby low-mass star forming regions. iButner et al.l (| 19951 ) found that X H i3co+ ranges from 



3 x 10" 11 to 4 x 10" 10 for dense cores with 2 x 10 21 cm 2 < Nn 2 < 1 x 10 22 cm 2 , where they 
estimated the H2 column densities from C 18 (J = 1 — 0) observations. For comparison, their values 
are indicated by crosses in Figure 01 Although their A H i3co+ tends to decrease with increasing 
A% 2 in the same way as for our data, their - Xm3m+ tend to be somewhat larger than the values 



predicted from eq. [6]. If -^h 13 co+ derived by lButner et al.1 (|1995) is fitted by a power-law with the 



same power index as in eq. [6], the coefficient becomes about twice that of eq. [6]. It is difficult, 
however, to judge whether this difference originates from the different environments because both 
the data have uncertainty at least by a factor of a few. If both the data are fitted by a power-law 
with a single index, then the best-fit function is given by 

logX H i3 CO + = (-10.89 ± 0.01) + (-0.732 ± 0.014) log(AH 2 /10 23 cm- 2 ) (7) 

with 1Z = 0.87. The mean fractional abundance (X H i3co+) is estimated to 1.72 x 10~ n with 
a standard deviation of 0.92 x 10~ n . In the present paper, we assume the constant fractional 
abundance of X^i3qq+ = 1.72 x 10~ n for the entire area. In future analysis, the spatial variation 
in fractional abundance will be taken into account. 



van Dishoeck et al.l (|1995l ) measured a smaller fractional abundance X H i3co+ ~ 0-75 x 10 
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toward IRAS 16293—2422, a young class bi nary YSO, located in the eastern streamer region of 



the p Ophiuchi molecular cloud complex (see lAndre et al.l 12000). They obtained the H2 column 



density of 2 x 10 23 cm~ 2 toward this YSO from the C 17 observations, somewhat larger than the 
mean H2 column density of p Oph evaluated from the 850 pm data. Their value, plotted in Fig. 0] 
by a filled square, is almost on our best-fit power-law of eq. [7]. Thus, we believe that our adopted 
fractional abundance of 1.72 x 10~ n is a plausible representative value for the p Ophiuchi main 
cloud. 

We note that several high-resolution interferometric observations of starless cores in relatively 
quiescent, low temperature star forming regions like Taurus have revealed that several gas-phase 
molecules including CO, HCO + , and their isotopes tend to deplete at the densities higher than 
10 5 cm 3 . One might interpret that the decrease in the H 13 CO + abundance with increasing column 
density is due to the effect of the depletion. However, the effect of the depletion should be small 
because of the following reasons. In cluster forming regions like p Oph, relatively high turbulent 
motions, protostellar outflows and stellar radiation tend to increase the temperature steadily and/or 
temporarily, preventing or slowing down the molecular depletion in starless cores. In addition, the 
high density part where the molecular depletion is significant is spatially localized and therefore 
the interferometric observations with high spatial resolution are needed to resolve the effect of 
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depletion in a cor e. However, our be am size (~ 20") is not enough to well resolve the structures 



inside each core (jAikawa et al.l 119951 ) . Furthermore, even if the depletion of gas-phase CO and 
HCO + is significant in p Oph, the fractional abundance of HCO + relative to CO does not change 
with increasing colu mn density because these carbon-bearing molecules tend to equ a lly de plete 
(Aik awa et al. 19951 ) . In c ontrast, the H 13 CO + fractional abundances of iButner et al.l ( 19951 ) and 



van Dishoeck et al 



(|1995l ). for both of which the H2 column densities are measured from the CO 
isotopes, do not show such a constant abundance, but decrease significantly with increasing column 
density, indicating that the dependence of the H 13 CO + abundance on the column density comes 
from the different physics from the molecular depletion. 



4.4. Distribution of Physical Quantities of Dense Cores and Comparison with Dust 

Cores 

In Table [IJ we present the physical properties of the 68 identified cores. In Table [21 we 
summarize the minimum, maximum, mean, and median values of each quantity. The histograms 
of the radius (i? C ore)) LTE mass (-Mlte)> an d mean density (n) of the H 13 CO + cores are shown in 
Figured) For comparison, the radius, mass, and mean density of the 850 /jm cores listed in Table 



7 of Jorgensen et al. (j2008l ) are indicated by the grey histograms in Figure [5j We selected only the 
850 /im cores located inside the same observed area as our H 13 CO + map. We also removed cores 
whose radii are smaller than the SCUBA beam size of 14" and corrected the radii of the other cores 



for the telescope beam size by using eq. [2] of llkeda et al.l (|2007l ). In this procedure, 43 dust cores 



were selected using the 2a noise level of 0.04 Jy beam -1 . The minimum, maximum, mean, and 
median values of each quantity of the 850 cores are summarized in Table [3j We note that as 
mentioned in in the 850 fim data the structures larger than 1.5' (~ 0.06 pc) are removed during 
the data reduction process and therefore the number of dust cores with radii larger than about 0.03 
pc is likely to be underestimated. 

The radius of the H 13 CO + cores ranges from 0.022 to 0.069pc and its distribution has a single 
peak at around the mean of 0.045pc (see Figure [5^). The maximum radius of the H 13 CO + cores 
is only about three times as large as the minimum radius. In other words, most of the cores have 
similar size. On the other hand, the distribution of the LTE mass is somewhat broad compared 
with that of the core radius. The LTE mass of the H 13 CO + cores ranges from 0.4 to 22 M Q and its 
distribution has a single peak at around the mean of 3.3M Q (see Figure [5b). The maximum mass 
is about 56 times as large as the minimum. The mean density of the H 13 CO+ cores also shows a 
similar broad distribution, ranging from 3.7 x 10 4 to 5.8 x 10 5 cm , with a single peak at around 
the mean density of 1.4 x 10 5 cm -3 (see Figure |5b). 

These distributions of the physical quantities of the H 13 CO + cores are qualitatively similar 
to those of the 850 /im cores, although the mean radius and mass of the H 13 CO + cores are about 
twice those of the 850 /jm cores. On the other hand, the mean density of the H 13 CO+ cores is 
smaller than that of the 850 ^m cores by a factor of 5. This large difference in the mean density 
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probably comes from th e fact that the 85 /im emission tends to trace the higher density parts of 
the H 13 CO + cores (see llkeda et al.l l2009). Another reason is due to the overlap effect. It is likely 
that dense cores often overlap one another along the line of sight. Although such overlapped cores 
can be separated in some degree for the H 13 CO + cores because of their velocity information, such 
separation is impossible for the two-dimensional clumpfind method that was used for the 850 fim 
map. In fact, the distribution of the core radius for the 850 pm cores has a tail toward the larger 
radius. A more prominent tail can be seen in the distribution of the core mass shown in Figure 
[5b, where the mean core mass is about 2.6 times as large as the median. The tail in the core 
mass, however, can be caused by the te mperature differe n ce. F or all the dust cores, the uniform 
temperature of T = 15 K is assumed in Ijprgensen et al.l (|2008l ). For several massive dust cores, 
however, the temperatures may be underestimated because they are located in the Oph A region 
where the temperatures are likely to be higher because of active star formation. 

The distribution of the FWHM line width of the H 13 CO + cores are shown in Figure El The 



line w idth ranges from 0.19 to 0.77 km s 1 with a mean of 0.49 km s . Following iMyers et al 



(|199ll ). we classified our cores into two groups: "thermal core" and "turbulent core", on the basis 
of the critical line width, dv CT , expressed by 



dv c . 



8ln2k B T 



1 



1 



1/2 



(8) 



^m H m ohS/ 

The critical line width of H 13 CO + is estimated to be 0.50 and 0.62 km s _1 for T = 12 K and 18 K, 
respectively. Figure [6] indicates that about 40 % of the cores are classified as the turbulent core. 
However, almost all the turbulent cores have transonic turbulent motions, an d highly turbulent 
cores (dvcore ^ 1 km s _1 ) found in massive star forming regions such as Orion A ( Ikeda et al.ll2007l ) 
are not seen in the p Ophiuchi main cloud. 



4.5. Comparison Between Cores with and without YSOs 



Recent observat ions of the Ophiuchi molecular cloud complex with the Spitzer Space Telescope 
(jPadgett et alj|2008l ) allow us to compare t he physical properties of H 13 CO + cores associated with 
and without YSOs (jJorgensen et al.1 120081 ). Previous studies have suggested that the physical 
prop erties of dense cores as s ociated with YSOs are different fro m those of cores without YSOs 
(e.g., iBenson fc Mverslll989l : Ijiiina et al.lll999l: IWalsh et alJ | 2007^ . Based on the spectral indices 



at the near- and mid-infrared wavelengths, iJorgensen et al.1 (120081 ) classified YSOs in L1688 into 4 
classes that represent the following evolutionary stages of YSOs: Class I, Flat Spectrum, Class II, 
and Class III. We use their YSO catalog to identify the H 13 CO+ cores with YSOs. We note that in 
their classification, the Class YSOs, the youngest objects, are included in Class I. For example, 
the prototypical Class object, VLA 1623, is classified as Class I in their list. 

The spatial distribution of the YSOs is shown in Figure which indicates that the distribution 
of the Class I and Flat Spectrum sources follows that of the H 13 CO + (J = 1 — 0) emission. On 
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the other hand, more evolved YSOs, the Class II and III sources, tend to be distributed in the 
entire observed area with no strong correlation with the distribution of the H 13 CO + ( J = 1 — 0) 
emission. Therefore, we expect that the Class I and Flat Spectrum sources are probably embedded 
in the dense gas, and use the two classes of YSOs to divide the H 13 CO cores into the following two 
groups: protostellar and starless cores. If either a Class I or Flat Spectrum object is located within 
the extent of a H 13 CO+ core, we regard the core as a protostellar core. The others are classified 
as starless cores, even if Class II or III sources are located within the cores. In this classification, a 
Class I or Flat Spectrum source is sometimes located within the extents of two or more cores which 
partly overlap one another along the line of sight. We classified all such cores as protostellar cores 
because we have no information on the line of sight velocities of the YSOs. Therefore, the number 
of protostellar cores is likely to be overestimated. According to the above procedure, we found 34 
protostellar and 34 starless cores. 

Figure [8] compares the physical properties of the starless and protostellar cores. The minimum, 
maximum, mean, and median values are summarized in Tabled! The radii (Figures EtO and masses 
(FiguresEb) of the protostellar cores tend to be larger than those of the starless cores. In particular, 
the mean LTE mass of the protostellar cores is about twice as large a s that of the starle ss cores. 



These results are consistent with those of the N2H + cores in NGC1333 (jWalsh et al.ll2007l ) and the 



H 13 CO+ cores in Orion A (jlkeda et al.l 120071 ). 



The line widths of the protostellar cores tend to be larger than th ose of the starles s core s 
(Fig uresEb). A simila r tendency is seen in the NH 3 cores observed by IBenson fc Myers 



and iJijina et al.1 (| 19991 ) , who interpreted that stellar winds and protostellar outflows create local 
turbulence. If stellar winds and protostellar outflows enhance internal turbulent motions for the 
protostellar cores, such cores must have had virial ratios smaller than the current values in the 
prestellar phase. Although the distribution of the virial ratios for the protostellar cores appears to 
be broader than that of the starless cores, there seems to be no significant difference in the virial 
ratio between the protostellar and starless cores (Figures Eti) : for both the starless and protostellar 
cores, the virial ratio has a mean of about 2 — 3, ranging from 0.4 to 8. Here, the dimensionless 
parameter a in eq.[2] is adopted to be unity, the value of a uniform sphere, for both the starless and 
protostellar cores. If for protostellar cores, the value of a centrally condensed density distribution, 
a = 5/3, is adopted, the virial ratios of the protostellar cores tend to be somewhat smaller than 
those of the starless cores that are expected to have flattened density distributions. This appears 
to contradict the idea that outflows enhance the virial ratios of the parent cores. About 60 % of the 
cores have virial ratios smaller than 2, for which the net internal energies are negative. Therefore, 
we conclude that the majority of the cores are likely to be self-gravitating. 

Figure Ek shows that the distributions of the mean densities of the protostellar and starless 
cores are similar to each other, although some of the protostellar cores have very large densities. 
Figure EF indicates that the starless cores tend to have larger aspect ratios than those of the 
protostellar cores. For example, for the starless cores, 5 (15%) out of 34 have aspect ratios larger 
than 2, whereas for the protostellar cores, no cores have such large aspect ratios. The starless cores 
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with aspect ratios larger than 2 have large virial ratios (a v ; r > 3) and are likely to be gravitationally 
unbound. 



4.6. Correlation among Physical Properties of Dense Cores 

4-6.1. Line Width-Radius and Mass-Radius Relations 



The line width measured from molecular line emission provides us with informa tion on kinetic 
properties such as thermal motion and turbulence. Based on data in the literature, iLarsonl (|198ll ) 
found that the line width of molecular clouds and cores is well correlated with the radius and 
mass, and the correlations are approximately of power-law form. He suggested that the power- 
law relations stem from processes of interstellar turbulence cascade. The strong correlations of the 
power-law relations also suggest that the molecular clouds and cores are nearly in virial equilibrium. 
Since his pioneering work, many studies have been carried out to investigate the dependence of the 
line width on its radius for various molecular clouds and cores using different molecular emission 
lines (e.g., ISanders et al.lll985l ; lDame et al.lll986l ; lHever fc Bruntll2004l ) . Observations based mainly 
on low-density trac ers such as 12 CO (J = 1 — 0) and 13 CO ( J = 1 — 0) have generally supported 



Larson's idea (e.g., iMyers 



1983 



Fuller Myersl ll992). In contrast, observations based on higher- 



density tracers such as C 18 0, CS, and H 13 CO + ha ve often shown that the dependence of the line 



width on its rad ius tends to be very weak (e.g., iTachihara et al. 



2002 



Ikeda et al. 



2007 



2009 



Saito et al 



200 



For massive star forming regions . the shallower power-law relations have also 



been found (|Caselli Mverdll995l : IPlume et al.lll997h . 



The line width-radius relation for the identified cores is presented in Figure [9^,. The filled 
squares and open circles indicate the starless and protostellar cores, respectively. As mentioned in 
§4.5( the protostellar cores tend to have larger line widths and larger radii than those of the starless 
cores. Such a tendency can be recognized in Figure [9^. The best-fit power-law for all the cores is 
given by 

log(cftw/km s" 1 ) = (0.319 ± 0.179) + (0.478 ± 0.131) log(i? corc /pc) , (9) 



with 7 Z = 0.410. As representat ive examples, the line width-radius relations derived by lLarson 



(|198ll ) and IHever &: Bruntl (120041) a re also plotted in Figure [9^ with the dotted and dashed lines, 
respectively. IHever &: Bruntl (|2004l ) applied the Principal Component Analysis to derive the line 
width-radius relation inside molecular clouds from 12 CO observations. The line widths of the 
H 13 CO + cores tend to be significantly larger than those of the Heyer & Brunt relation and somewhat 
smaller than those of the Larson relation. We note that the line widths of the Larson relation tend 
to be larger over t he scale 0.01 — 0.1 pc compared to the line width-radius relations derived by 



other authors (e.g.. ISolomon et al. 



1987 



Myerall983l ) and thus may correspond to the upper b ound 



of the relations . Altho ugh the power index of our cores is not far from those of ILarsonl (Il98ll ) and 
Hever &; Bruntl (|2004l ). the correlation is ver y weak. Such a weak correlation is also seen for the 
H 13 CO + cores in Orion A (jlkeda et al J 120071 ). In £[5]we will discuss the line width-radius relation 
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of the H 13 CO + cores in more detail. 



To separate the turbulent motion from the thermal motion within each core, we plot in Figure 
[9b the nonthermal line width (1v-nt[= (dv% ot — 8 In 2 ksT / irniH) 1 ! 2 } against core radius. The best-fit 
power-law for all the cores, shown by the solid line in Figure [9b, is given by 

log(dt; NT /km s^ 1 ) = (0.501 ± 0.218) + (0.634 ± 0.159) log(i*/pc) (10) 

with TZ = 0. 440. For compa r ison, the nonthermal line width-radius relation for low- mass cores 
compiled by iCaselli fc Mversl ()1995l ). (c^NT/km s" 1 ) = 1.2(i t!/pc) 0,53 , is indicated by the dashed 
line, where we replaced the radius of ICaselli &, Myers! ()1995l ) (Rcm) by the radius of our defini- 
tion. See the Appendi x in detail. Although our relation has a power-law index similar to that of 



Caselli &; Mversl (119951 ). the coefficient is somewhat larger. This suggests that turbulence within 



our cores is not as much dissipated as the Caselli & Myers relation predicts. This tendency is 
different from that of the H 13 CO + cores in Ori on A where the dv^T-R COTe relation agrees well with 
the Caselli & Myers relation (jlkeda et al.ll2007l ). In £}5] we will compare our cores with the Orion A 
cores and discuss this point in more detail. 



Recently, IVolgenau et all (120061 ) attempted to measure the relationship between the line width 
and size inside three protostellar cores having multiple pro tostars in the Perseus molecular cloud, 
using the Regions-Of-Contrast (ROC) analysis proposed by lOstriker et all (|200ll ). Our line width- 
radius relation is essentially different from that derived from the ROC analysis, and therefore, it 
is difficult to directly compare the two relations. In the ROC analysis, the size is just a region 
size that is not related to a dense core. When the ROC analyis is applied to substructures inside 
a core, the line width-size relation tends to ha ye a large scatter, esp ecially when the objects have 
complicated structures like those observed by IVolgenau et al.l (|2006l ). in which the disk rotation, 
disk instability, and stellar feedback such as stellar winds and outflows influence the local structures 
significantly. In contrast, in our analysis, the size (or radius) corresponds to that of an identified 
core. Therefore, the relatively large scatter in our line width-radius relation is likely to reflect the 
difference in the dynamical states of the cores (see also §5.2h . In fact, recent numerical simulations 
of turbulent molecular clouds demonstrated that the cores identified directly from the 3D density 
distribut ion are in various dynamical states and s how similar large scatter in the line width-radius 
relation (|Nakamura fc Lilhoosl : foffner et alJboQsl ). 



In Figure [9b, we plot the LTE mass against core radius. The best-fit power-law is given by 

log(M LTE /M ) = (3.95 ± 0.34) + (2.64 ± 0.25) log(R coie /pc) (11) 

where TZ = 0.790. For comparison, the relation of Mlte = 47r/i7Bffn.R 3 ore /3 is indicated by the 
dashed line in Figure [9b, where n = 1.35 x 10 5 cm~ 3 . The power-law index of ~ 3 indicates 
that the mean density of each core is comparable to the critical density of H 13 CO + (J = 1 — 0), 
n CT ~ 8 x 10 cm" 3 . T he coefficient of our MLTE--Rcore relation is an order of magnitude larger than 
that of iLarsonl ( 198ll ). which is mainly based on 13 CO (J 
10 3 cm 



2009|). 



1 — 0) whose critic al density is about 
3 . This tendency is in good agreement with that of the Orion A cores ( Ikeda et al J 120071 . 
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4-6.2. Virial Ratio 



The boundedness of a core is often estimated using the virial ratio. In Figure [TU1 we plot against 
LTE mass the virial ratio calculated using eq.[2], where for both the protostellar and starless cores 
a dimensionless parameter a is set to unity, the value of a uniform sphere. The best-fit power-law 
of our cores is given by 



loga vir = (0.528 ± 0.018) + (-0.600 ± 0.057) log(M LTE /M ( 



(12) 



with 1Z = 0.817. The virial ratio tends to increase with decreasing the LTE mass, and has a mean 
of 2.4, ranging from 0.4 to 8. More than a half the cores have virial rati os smaller than 2 and thu s 
the majority of the cores are more or less self-gravitating. According to lBertoldi &; McKed ([1992), 



a self-gravitating core confined by ambient pressure ha s a virial ratio of a v \ r = 2. 06(Mlte/-Mj) 2 / 3 



where Mj is the Jeans mass defined by eq. [2.13] of iBertoldi &: McKed (|1992l ). For comparison, 
the virial ratio of a self-gravitating core confined by ambient pressure is plotted by the dashed line 
in Figure fTUl Here, using the velocity dispersion of 0.35 km s _1 and the ambient gas density of 
10 5 cm -3 (see the derivation of these values in $5j), the Jeans mass is estimated to be about 2.4 
M©, almost equal to the median mass of the identified cores. The best-fit power-law is very close 
to the virial ratio of a self-gravitating core confined by ambient pressure. The virial ratios of the 
identified cores tend to be, however, smaller than that of a non self-gravitating, pres sure-confined 
core with a v ; r = 2.9(Mlte/-^j) _2/ ' 3 ! for a given core mass ( Bertoldi McKed Il992l ). Therefore, 
for our cores, both the self-gravity and ambient pressure play an important role in dynamics of the 
cores. In $5]we will discuss the dynamical states of our cores in more detail. 



Discussion 



5.1. Core Mass Spectrum in the p Ophiuchi Main Cloud 



The mass spectrum of our cores is plotted in Figure [TTJ The core mass spectrum appears to 
be fitted by a two-component power-law. There seems to be a break at around 7M@, about twice 
the mean LTE mass of 3.4M Q . Above the break, the mass spectrum is steeper, while below the 
break, it is flattened. The mass spectrum can be fitted by the following two-component power-law: 



dN/dM oc M-°- 4d±0 - 03 for M LTE < 7M 
« M -2.4±o.3 for Mlte > 7M& _ 



(13) 
(14) 



This mass spectrum is broadly consistent with those of the dust cores in p Oph, alt hough for the dust 



cores the low-mass parts are somewhat steeper and the break masses are smaller (IMotte et al 



1998 



Johnstone et alJboOOa : Stanke et all2006 ; Stamatellos et all2007 ; Enoch et alJbood : Simpson et al 
2008). For example, Stamatellos et al. ( 20071 ) revised the mass spectrum of the dust cores identified 



by lMotte et al.l (119981 ) by reestimating the dust temperatures. The revised mass spectrum is fitted 
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by a two-compon ent power-law with a low-mass index of —1.5, a high-mass index of —2, and a break 
at around 1M@. ISimpson et al.l (j2008l ) reanalyzed the SCUBA 850 pm archive data and obtained 
the core mass spectrum that can be fitted by a three-component power-law with a low-mass index of 
—0.7, intermediate-mass index of —1.3, high-mass index of —2.35, and two breaks at around O.7M 
and 2Mq. We note that the slopes and the brea k masses of the co re mass spectra are affected by 



the number of bins, if it is not appropriately set (jRosolowskyl 120051 ) . Thus, the cumulative form of 
the core mass spectra is sometimes used when th e number of cores is small, although it requires 
a more complicated uncertainty analysis (see e.g., iReid fc Wilson! 120061 ) . As for the parameters of 
the clumpfind, the slopes are almo st independent of the thr eshold and stepsize of the clumpfind as 
long as they are appropriately set (jlkeda fc Kitamurall2009l ). For the core mass spectra mentioned 
above, the difference in the break mass is likely to come from the fact that the H 13 CO+ emission 
traces more extended, less dense structures than the dust continuum emissions. The break mass 
may also be affected by the "confusion" which means the situation that multiple cores with similar 
vlsr are overlapped along the same line of sight and t herefore cannot be separated into individual 
cores using the clumpfind (see llkeda et alJl2007l . 12009 ) . 



Based on the ISOCAM observations, iBontemps et al.l (|200ll ) identified over 200 YSOs in p 
Oph and derived the mass function of 123 Class II YSOs that is well fitted by a two-component 
power-law with a low-mass index of —1.15, a high-mass index of —2.7, and a break at around 
0.55A/q. The mass spectrum of the H 13 CO + cores is roughly similar in shape to the stellar IMF in 
p Oph, although the slopes in the low mass and high mass parts of the H 13 CO + core mass spectrum 
are somewhat shallower and the break mass is one order of magnitude larger. 



5.2. Dynamical State of H 13 CO + Cores in the p Ophiuchi Main Cloud 



Virial theorem is useful for analyzing the dynamical states of dense cores. The virial equation 
for a uniform spherical core is given by 

Id 2 / 

-— = U + W + S (15) 

where the terms, I, U, W, and S, denote the moment of inertia, internal kinetic energy, gravitational 
energy i ncluding the effect of magnetic field, and surface pressure, respectively, and are given as 
follows jNakanolll998h : 



U 
W 



3Mdvl t 
81n2 
3 GM 2 



5 RcOTC 



(16) 

(17) 
(18) 



The values $ and <3? cr are, respectively, the magnetic flux penetrating the core and the critical 
magnetic flux above which the magnetic field can support the core against the self-gravity. Here, 
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the core is assumed to have condensed from much lower density medium and the radius of the 
magnetic flux tube penetrating the core is much smaller than that before contraction [see §2 of 
Nakanol 41223)]- 

P ex is the surface pressure including both thermal and turbulent components. M 
is the core mass and is chosen to be the LTE mass of the core. 

All the above terms in the virial equation except P ex and $ can be estimated from the physical 
quantities listed in Table [T] It is difficult to estimate the surface pressures exerted on individual 
cores directly from our data. Instead of deriving the surface pressures of individual cores, we adopt 
a representative surface pressure from the average densities and velocity dispersions that were 
measured in subclupms Oph A, B, C, E, and F: (P ex ) ~ (p) (c) 2 - Here, (p) is the average density 
and (ex) is the average velocity dispersion including both the thermal and turbulent components as 
(a) 2 = (cnt) 2 + c 2 ., where (<tnt) and c s are the velocity dispersion of the nonthermal component 
and the isothermal sound speed, respectively. The average density and velocity dispersion were 
evaluated as (p) /(2.33m//) « (0.5 — 1) X 10 5 cm -3 and (a) « 0.3 — 0.4 km s _1 , respectively. The 
average surface pressure is then given by (Pex) /&B ~ 3 x 10 6 K cm -3 . This value is about a few 
times as large as the thermal pressure [~ 10 5 cm" 3 x (12 — 18) K ~ (1 — 2) x 10 6 K cm -3 ]. In the 
following, we adopt the above (P ex ) as a representative value for our cores. This surface pressure 
is about a half the average internal pressure inside the cores («6x 10 6 K cm -3 ) and close to the 



lower limit of the critical ambient pressure for the dust cores obtained by iJohnstone et al.l (|2000al ) . 



On the other hand, there is only one reliable mea surement of the magnetic field strength toward 
the p Ophiuchi main cloud. ICrutcher et al.1 (119931 ) performed OH Zeeman effect measurements 
toward two positions in p Oph with a beam size 18'. For one of the two positions that well covers 
the entire observed area of our data, they derived the line-of-sight magnetic fi eld strength of abou t 
10 pG at the low density cl oud envelope of Nr 2 « 5 x 10 21 cm -3 (see also iTroland et al.lll996l ). 
Based on this measurement. ICrutcherl (|1999l ) estimated the magnetic flux normalized to the critical 
value to be ^/^cr — 0.4. If we adopt this value as a representative value of $/$ CT for our cores, 
the magnetic effect is likely to be minor: it reduces the gravitational energy term only by 16%. 
However, recent turbulent simulations have demonstrated that cores formed out of turbulent clouds 



can have much larger values of ^/^cr than the cloud initial values (e.g. JDib et al.ll2007l ) and thus it 
is very difficult to assess the values of $/$ C r for the individual cores without direct measurements 
of the magnetic fields associated with the cores. In the following, we simply assume &/<& C r = 0. 

The equilibrium line, U + W + S = 0, is shown by the solid line in Fig. [T2l where the surface 
term (S) is plotted against the gravitational energy term (W), both normalized to the internal 
energy term (U). For the cores that lie below the solid line, the value of U + W + S is negative 
and thus expected to be bound. All the others are unbound and expected to disperse away, if 
they do not gain more mass through accretion and/or merging with other cores, or reduce internal 
support through turbulence dissipation. More than a half the cores lie below the equilibrium line, 
and are thus bound: they are expected to collapse and form stars. This is also true even in the 
presence of magnetic field as long as the magnetic flux does not exceed about a half the critical 
value. Furthermore, the majority of the cores lie below the line of S = W (dashed line), indicating 



-20- 



that the surface term is more important than the gravitational energy term. Even for the cores 
lying above the line S = W, the surface pressure appears to be dynamically important for almost 
all such cores because the deviation from the line is small. We note that this plot should be valid 
statistically, but for any individual points the true surface pressures could be higher or lower than 
the representative surface pressure. 

Our virial analysis indicates that the formation and evolution of dense cores in the cluster 
environment are likely to be strongly influenced by the external compression due to local turbulent 
motions (see also Dobashi et al. 2001 for importance of external pressures at cloud scales). Such 
local compression may be responsible for the formation of binary and multiple stars or formation 
of substellar objects (brown dwarfs or planetary mass objects) in each core under the cluster 



environment (Hennebelle et al 



2003 



Gomez et al 



2007 



Whitworth et al.1 120071 ) . 



5.3. Comparison with Orion A 

Th e Orion A molecula r cloud is the nearest giant molecular cloud located at a distance of 
480 dc iomzdet al.1 fl98lh. about 4 times as distant as the p Ophiuchi main cloud. Recently, 



Ikeda et al.1 (|2007h carried out the H 13 CO + ( J = 1 — 0) core survey in the whole region of Orion 
A using the Nobeyama 45 m telescope and identified 236 cores by the same method as adopted 
in the present paper. Their core sample is ideal for comparison with our cores because the same 
molecular emission line, the same telescope, and the same core identification procedure are used. 
In this subsection, we compare the physical properties of our cores in p Oph with those in Orion 
A. 

Figures [TBI and [T4l compare the physical properties of the p Oph cores with those of the Orion 
A cores. The open squares and crosses indicate the p Oph cores and the Orion A cores, respectively. 
The mean radius of the p Oph cores, 0.045 pc, is about three times smaller than that of the Orion 
A cores, 0.14 pc. This probably reflects the different spatial resolutions (or different distances) for 
both the observations. In fact, the area where the p Oph cores are distributed is well separated 
from that of the Orion A cores on the plots of the line width-radius and mass-radius relations (see 
Figure [TBI) . On the other hand, the mean mass of the p Oph cores, 3.35 M & , is only about 3.6 
times smaller than that of the Orion A cores, 12 Mq. This is because the mean density of the 
Orion A cores, 1.6 x 10 4 cm -3 , is much smaller than that of the p Oph cores, 1.4 x 10 5 cm -3 . In 
particular, the mean density of the Orion A cores is smaller than the critical density of the H 13 CO + 
(J = 1 — 0),8xl0 4 cm -3 , and therefore the Orion A cores are likely to contain substructures that 
could not be spatially resolved. 

The mean line width of the p Oph cores, 0.488 km s" 1 , is almost the same as that of the Orion 
A cores, 0.52 km s _1 , in spite of the different core radii. The ranges of the line widths are also 
similar to each other (see Figure [TB]) . On the other hand, the virial ratio of the Orion A cores seems 
three times larger than that of the p Oph cores for a given LTE mass, as shown in Figure HH This 
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difference may be explained also by the different spatial resolutions as follows. The virial ratio is 
proportional to both the radius and the square of the line width for a given LTE mass. Since the 
line widths of both the cores are comparable to each other and the radii of the Orion A cores are 
roughly three times larger than those of the p Oph cores, the virial ratio of Orion A becomes three 
times larger than that of p Oph for a given LTE mass. 

To elucidate whether or not the different spatial resolutions dominate the different core prop- 
erties in the two clouds, we try to eliminate the effect of the different spatial resolutions as follows. 
We smoothed the original p Oph data on a coarser grid with a grid spacing of 80" by convolving 
with a 2D Gaussian kernel with a FWHM of 80" in the a-5 space. The smoothed data correspond 
to the data that would be observed by the Nobeyama 45 m telescope if the p Ophiuchi main cloud 
were located at the same distance as the Orion A molecular cloud. Then, applying the clumpfind 
to the smoothed data, we reidentified 16 dense cores, which are plotted by the filled circles in 
Figures [TBI and [T4l In Figures [T3"a and[T3b. the p Oph cores identified from the smoothed data are 
located just in the areas where the Orion A cores are distributed. In contrast, the virial ratios of 
the reidentified p Oph cores are somewhat smaller than those of Orion A for a given LTE mass (see 
Figure [Lib'). T h e larg er virial ratios for the Orion A cores may reflect the larger Xr13qq+ adopted 
bv llkeda et al.1 ([20071 ). 4.8 x 10 -11 , which is derived using the 13 CO abundance measured in the 
cloud envelope of Tauru s, which has much low er column densities than Orion A. If we use the 13 CO 



abundance obtained by iFrerking et alJ 



1982) for p Oph, which has column densities comparable 



to the region observed bv llkeda et al.l ([20071 ). the H 13 CO + abundance is reduced by a f actor of 



2, close to our value. Note that our value is in agreement with that of llkeda et al.1 (120071 ) within 
uncertainty mentioned by them. Using our value of X H i3co+i the distribution of the Orion A cores 
in the virial ratio-LTE mass plot becomes consistent with that of the p Oph cores identified from 
the smoothed data, as presented in Figure [T4"b. Although in that case, the LTE masses of the Orion 
A cores shift upward in Figure [13b, the distributions of the p Oph and Orion A cores in the LTE 
mass-radius relation plot are still well overlapped with each other. Therefore, we conclude that 
there are no clear differences in the core properties between the p Oph and Orion A clouds. The 
apparent differences in the core properties between the two clouds are likely to be caused mainly 
by the different spatial resolutions. 



5.4. Line Widths of the H 13 CO + cores and Implication for Turbulent Generation 



Figure [TBI shows that there are no highly turbulent cores having dv coie > 1 km s" 1 in p Oph. 
According to Ikeda et al. (j2007l ). all the highly turbulent cores in Orion A are located within 1 
pc from the M42 HII region and are probably influenced greatly by the nearby OB stars. The 
tu rbulent cores are expected to be responsible for massive star formation (see the discussion of §4.3 
of llkeda et al.ll2007l ). The absence of such highly turbulent cores in p Oph may imply that massive 
O stars will not form in p Oph under the current environment. 



It is also interesting that the range of the line widths of the p Oph cores identified from the 
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smoothed data almost coincide with that of the original p Oph cores. If we fit a single power-law 
to both the p Oph cores identified from the original and smoothed data, the best-fit result becomes 
dv core = (0.599±0.104)i?[! o l Jg 3= ' =0 ' 060 with a very small correlation coefficient of 0.137, which is shown 
by the solid line in Figure [T3j In other words, the line width appears to be almost independent of 
the core radius. Since the cores identified from the smoothed data contain adjacent several cores 
identified from the original data, the nearly independent relationship between the line width and 
core radius, or the nearly flat line width-radius relation, suggests that the inter-core motions among 
the neighboring cores are al most compar a ble to the internal motions in the individual cores. Such 
a feature is pointed out by I Andre et al.l (|2007l ) who measured the velocity difference among the 
neighboring cores from the centroid velocities of the cores observed by N2H + ( J = 1 — 0). If the line 
width-radius relation is flat (dv corc « const.) and the core mass is proportional to -R 3 ore , then the 
virial ratio is scaled as a v \ r oc R cme dv1 ore l '(GMlte) oc M l ^Z 3 , a similar power-law to that derived 



bv lBertoldi fc McKed (|1992| ). This again suggests the importance of ambient turbulent pressure in 
dynamics of the cores as discussed in §5.21 



The line width-radius relation of the p Oph cores appears to be inconsist ent with one of t h e mos t 
reliable measurements of the line width-size relation, recently obtained bv lHever &: Bruntl (|2004l ). 
who found a strong correlation between the line width and size based on the 12 CO observations 
toward 27 nearby giant molecular clouds, despite the large differences in cloud environments and 
local star formation activity. Their relation, plotted by the dashed line in Figure [T3l [(d^/kms -1 ) = 
l.l(i?/pc)° ], leads to the interpretation that most of the tur bulent energy comes fr om the largest 
scale that is comparable to the cloud size (see the Appendix). iHeyer Bruntl (|2007l ) claimed that 
the driving mechanisms of turbulence at small scales such as protostellar outflows and stellar winds 
may not play an important role in dynamics of molecular clouds. The FWHM line width (AV ~ 1.5 
km s _1 ) and radius (R ~ 1.5 pc) measured from the 13 CO (J = 1 — 0) line toward the whole p 
Ophiuchi main cloud are in good agreement with the Heyer & Brunt relation. In contrast, the p Oph 
cores identified from the original data tend to deviate upwards from the Heyer & Brunt relation. 
The nearly flat line width-radius relation of the p Oph cores may also suggest that turbulent energy 
has been injected at the scales smaller than R ~ 0.35 pc, at which our best-fit power-law intersects 
with the Heyer & Brunt relation. This scale is in reasonable agreement with the chara cteristic length 
scale of the outflow-d riven turbulence estimated from the theoreti cal consideration (IMatznerl 120071 ; 
Nakamura fcLil20n7h . Based on 13 CO and C 18 observations, ISwift fc Welch! feooah derived a 
similar characteristic scale of the outflow-driven turbulence of ~ 0.05 pc toward L1551. However, 
the universality of the Heyer & Brunt relation implies that the outflow-driven turbulence may be 
limited to the localized dense regions where active star formation occurs. A caveat is that it is 
unclear whether the Heyer Sz Brunt relation, based on the 12 CO observations, is still applicable 
to dense gas where stars are forming. Therefore, the deviation of the line widths for p Oph from 
the Heyer-Brunt relation may be apparent. However, it is worth noting that the nearly flat line 
width-radius relation is consistent with the virial ratio-mass relation of the pressure-confined, self- 
gravitating cores (a v i r oc A/- 2 / 3 ). This may suggest that the characteristic length scale of 0.35 pc 
is related to the formation of dense self-gravitating structures. In any case, if such a small scale 
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driving of turbulence is common in the cluster environment, star formation in clusters is likely to 
continue over several dynamical times. 

On the other hand, the line widths of our cores deviate downward from the Larson relation 
which may correspond to the upper bound of the observed line width-radius relations. Until now, 
many authors have derived the line width-radius relations for various molecular clouds and cores. 
Those studies show a scatter in the index and coefficient of the power- law line width-radius relations. 
The line widths of the Larson relation, (Au/kms -1 ) = 1.9(i?/pc) ' 38 , tend to be somewhat larger 
over the range of 0.01 — 0.1 pc compared to the other relations and thus may correspond to the 
upper bound of the relations. For exampl e, on the basis of the Massachusetts-Stony Brook Galactic 



Plane Survey data, ISolomon et al.l (j 19871 ) derived the relation of (Av/kms ) = 2.1(R/pc) , no t 



far from the Heyer & Brunt relation. In nearby low-mass star forming regions iFuller &; Myersl (jl992l ) 
compiled the data of dense cores and derived the power-law relation of (Au/kms -1 ) = 1.0(i?/pc)° 
for starless cores. If the line width-radius relation of the low-density gas in p Oph follows the Larson 
relation, instead of the Heyer & Brunt relation, the smaller line widths of the H 13 CO + cores may 
suggest that the cores have formed preferentially in regions where supersonic turbulence dissipated 
and thus the cores so formed have smaller line widths compared to the Larson relation. If this is 
the case, the nearly independent relationship between the line width and core radius implies that 
the turbulent field at the small scales has not been relaxed yet, and therefore the timescale of core 
formation may be of the order of a dynamical time. In that case, it is unclear how protostellar 
outflows contribute to maintain supersonic turbulence in the cluster environment. 



6. Summary 

Using the archive data of the H 13 CO + (J = 1 — 0) molecular line emission taken with the 
Nobeyama 45 m radio telescope, we analyzed the molecular gas distribution in the central dense 
region of the p Ophiuchi main cloud. We summarize the primary results of the paper as follows: 

1. We compared the global distribution of the H 13 CO + emission with that of the 850 pm 
dust continuum emission and found that the overall spatial distributions are similar to each other, 
while the H 13 CO + emission appears to cover a larger area than the dust continuum emission. By 
comparing between the H 13 CO+ and 850 pm maps, we revealed that the fractional abundance of 
H CO + relative to H2 decreases with increasing the H2 colum n density as X w i3nn + oc N n2 ' 



This tendency is consistent with the theoretical prediction (e.g., iNakano et al.ll2002l ). The mean 
fractional abundance in the region is estimated to be 1.72 x 10 . 

2. From the 3D data cube of the H 13 CO+ emission we identified 68 dense cores using the 
clumpfind method. From comparison with the positions of YSOs recently identified by the Spitzer 
space telescope, the cores are classified into the following two groups: 34 protostellar and 34 starless 
cores. The radii, masses, and line widths of the protostellar cores tend to be larger than those of the 
starless cores. The virial ratio tends to increase with decreasing the LTE mass, although we found 
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no significant difference in the viriai ratio between the protostellar and starless cores. Furthermore, 
the viriai ratios of the cores can be well described by the model of a self-gravitating core confined 



by ambient pressure derived by iBertoldi McKed (|1992l ). suggesting that both self-gravity and 



ambient pressure play an important role in dynamics of the cores. 

3. The mass spectrum of the H 13 CO + cores can be fitted by a two-component power-law that 
resembles the stellar IMF in p Oph. The core mass spectrum has a break at around 7M & , close to 
the mean LTE mass of the cores. It implies that the stellar IMF may be at least partly determined 
by the core mass distribution. 

4. Applying the viriai analysis, we conclude that most of the cores are bound and are expected 
to collapse. Furthermore, for the majority of the cores, the surface pressure term is more impor- 
tant than the gravitational energy term. This result suggests that in the cluster environment the 
formation and evolution of the cores may be regulated largely by the surface pressures. 

5. We compared the physical properties of the H 13 CO + cores in p Oph with those in Orion 
A. To eliminate the effect of the different distances, we smoothed the original p Oph data into a 
coarser grid so that the spatial resolutions coincide with each other. We identified 16 cores from 
the smoothed data using the clumpfind. These cores would be observed by the Nobeyama 45 m 
telescope if the p Ophiuchi main cloud were located at the same distance as the Orion A molecular 
cloud. The physical properties of the p Oph cores identified from the smoothed data appear to 
resemble those of Orion A. The fact that the mean density of the Orion A cores (1.6 x 10 4 cm~ 3 ) 
is somewhat smaller than the critical density of H 13 CO + ( J = 1 — 0) transition of 8 x 10 4 cm~ 3 
suggests that the Orion A cores contain clumpy substructures that could not be spatially resolved. 

6. The range of the line widths of the p Oph cores identified from the smoothed data almost 
coincides with that of the original p Oph cores. This suggests that the line width-radius relation is 
more or less flat over the range of 0.01 — 0.1 pc for the dense gas detect ed by H 13 CO + . Such a flat 



relation is inconsistent with that derived from the 12 CO observations bv lHever fc Bruntl ( 2004 ) who 
found that the line width correlates strongly with the size for the low-density molecular gas. At the 
scales below 0.3 pc, the line widths of the p Oph cores deviate upwards from the Heyer & Brunt 
relation. This may be due to turbulence driven by protostellar outflows. However, if the actual line 
widths of the low-density gas in p Oph agree with the values of the Larson relation, then the nearly 
independent relationship between the line width and core radius may be interpreted as evidence 
that the cores have formed preferentially in regions where supersonic turbulence dissipated. 



A. Comparison with Other Line Width-Radius Relations 



In the present paper, the FWHM line width and core radius defined by our clumpfind method 
are used for showing the line width-radius relation of the iden tified cores. Since the different authors 



use the different definitions of the line width and core radius (He ver Brunt 



2004 



1995 



Larson 



1981 



Caselli Myers 



we here convert the other line width-radius relations using our definitions. 
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Hever Brunti (|2004j) applied the PCA analysis to derive the relationship between the line 
wdith and the size inside the molecular clouds. Their line width A^hb and size Lhb correspond to 
the autocorrelation widths on both the spatial and the velocity axes. In the following, we transfer 
them into our FWHM line width and radius to do a direct comparison with our result. 

If both a line and intensity profiles of a core can be ap proximated by Gaussian profiles, 
then the line width and size defined by iHever &: Brunti (|2004l ) are given by Av hb = 2<7„ and 
Lhb = 2a /- ; , where tr „ and ul are the dispersions of the line and intensity profiles, respectively (see 
Hever & Brunthooi l. On the other hand, the FWHM line width is defined as Awfwmh = V8 In 2a v . 
Thus, the FWMH line width can be rewritten as Ai>fwhm = 1-18Avhb- 

In the present paper, the core radius, R CO re, is determined by taking the projected area enclosed 
by the 2 a level contour and by computing the radius of the circle required to reproduce the area. 
For a core with a Gaussian profile, our core radius is given by i? CO re = ^hb \/ m C^peak T^th ) > where 
Tpeak an d T t h are the peak antenna temperature and the threshold temperature of the 2 a level, 
respectively. From Table HJ the mean peak antenna temperature of the 68 identified cores is 
estimated to be (T pca k) = 1.14K. If we adopt (T pea k) as a representative value, i? C ore = 0.91Lhb- 



According to IHever fc Brunti (120041 ). the linewidth-size relation for GMCs can be fitted by 
the power-law of (A-UHB/km s" 1 ) = 0.87(Lhb/pc) ' 65 . Using Aufwhm = 1.18Adhb and R core = 
0.91Lhbj the Heyer & Brunt relation can be rewritten as (AupwHM/km s _1 ) = 1.09(i? core /pc) ' 65 , 
which is shown by the dashed line in Figure [T3k . 



Caselli Sz Myers (| 19951 ) also used the different definition for the core radius. They used the 



half-maximum contour to define the core size. Assuming that a core has a Gaussian shape, the 
half-maximum radius is equal to -Rcm = y^n 2/ ln(T pea k/T t h)i? C ore) and Rqm = 0.65R core for the 
mean peak antenna temperature of (T pea k) = 1.14K. We thus replace the original Caselli & Myers 
relation, (dv NT /km s" 1 ) = 1.51(i? C M/pc) - 53 , with ((fo NT /km s -1 ) = 1.2(i? coro /pc) - 53 . 



Larson! ()198ll ) used the 3D velocity dispersion (030) and the maximum core size (L max ). As- 
suming that -L max = 2R cove , his original relation, (o"3D/km s _1 ) = l.l(L max /pc) ' 38 ), is rewritten 
as (Av FW HM/km s _1 ) = 1.9(i? core /pc) - 38 . 
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Fig. 1. — (a): H 13 CO + (J = 1—0) total integrated intensity map in the velocity range of "Ulsr = 1 — 6 
km s -1 toward the p Ophiuchi main cloud. The grey scale indicates the integrated intensity in units 
of K km s _1 . The contours start from 0.3 K km s , corresponding to the 3 a noise level, at intervals 
of 0.2 K km s _1 . (b): 850/mi image of the same area as in panel (a), obtained with the SCUBA 
at the JCMT. The grey scale indicates the intensity in linear scale from -0.4 to 4 Jy beam . 
The contours start from 0.2 Jy beam -1 at intervals of 0.2 Jy beam -1 . The dense subclumps are 
designated by A, Bl, B2, C, E, and F. 
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Fig. 2.— H 13 CO + (J = 1 - 0) total integrated intensity maps taken with the Nobeyama 45 m 
telescope toward (a) the Oph A region and (b) the Oph Bl, B2, and C regions. The grey scale and 
contours are the same as those of Figure [Ik. In panel (a), the two dahed lines indicate the position 
of the filaments seen in the 850//m map (jWilson et al.l ll999). The positions of the submillimeter 
sources, SMI, SM1N, and SM2, and the prototypical Class object, VLA 1623, are indicated by 
crosses. In each panel the arcs and holes discussed in the text are indicated by numbers. The beam 
size of the Nobeyama 45 m telescope is shown at the lower left corner of each panel. 



-33- 




Fig. 3. — Three-dimensional representation of the antenna temperature (T|) in the a-8-v^sK space 
toward (a) the Oph A region, (b) the Oph Bl, B2, and C regions, and (c) the whole observed area. 
For panels (a) and (b), the areas projected on the plane of the sky are the same as those of Figures 
[2ji and [2b, respectively. The color shows the iso-temperature surfaces. The color bar indicates the 
antenna temperature in units of K. In each panel the numbers have the same meaning as in Figure 
[2j (Mpeg animations of Figs. [2^, [2b, and [2k are available in the online journal.) 



10 



12 



10 



21 



10 



11 



10 



11 



10 



24 



N (cm" 2 ) 



Fig. 4. — Fractional abundance of H 13 CO + relative to H2 against the H2 column density in the p 
Ophiuchi main cloud. The solid line shows the best-fit power-law fo r the p Oph data. For compari 



son, t he fractional abundances of H 13 CO + derived by lButner et alJ (1993) and Ivan Dishoeck et al 



(|1995l ) are plotted by crosses and a filled square, respectively. The dashed line shows the best-fit 
power-law for both their data and ours. 
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Fig. 5. — Histograms of (a) the radius, (b) core mass, (c) mean density of the H 13 CO + cores (open 
histograms) and the 850 fim dust cores (grey histograms) in the p Ophiuchi main cloud. 
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Fig. 6. — Histogram of the FWHM line width of the H 13 CO + cores in the p Ophiuchi main cloud. 
The open and grey histogams indicate the thermal and turbulent cores, respectively. Note that 
the gas temperature is assumed to be 12 K for all the cores except for those located in the Oph A 
region, for which T = 18 K. 
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Fig. 7. — Distribution of YSOs identified with the Spitzer Space telescople. Contours indicate the 
total integrated intensity of the H 13 CO + (J = 1 — 0) emission, as shown in Figure [lb. The red and 
blue circles are for Class I sources including Class and Flat Spectrum sources, respectively. The 
green and cyan squares are for Class II and Class III sources, respectively. The cores associated 
with the circles are defined as protostellar cores in this study. 
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Fig. 8. — Histograms of (a) the radius, (b) LTE mass, (c) FWHM line width, (d) virial ratio, 
(e) mean density, and (f) aspect ratio of the H 13 CO + cores. For each panel, the open and grey 
histograms are for the cores without and with YSOs, respectively. 




Fig. 9. — (a) Line width-radius relation, (b) Nonthermal line width-radius relation, and (c) Mass- 
radius relation of the H 13 CO + cores. The filled squares and open circles are for the cores without 
and with YSOs, respectively. In each panel, the best-fit power-law is indicated by the solid line. 
In panel (a) the Larson and Heyer & Brunt relations are plotted by the dotted and dashed lines, 
respectively. Note that we extrapolated these relations down to 0.01 pc, for comparison. In panels 
(b) and (c), the dashed lines indicate the Casern &; Myers relation and the relation of Mlte = 
47r//m#nit!j? OI . c /3 (oc -Rc 0re ), respectively 
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Fig. 10. — Virial ratio-mass relation of the H 13 CO + cores. The filled squares and open circles 
are the same as those in Figure [9j The virial ratio of a self-grav itating core confined by ambient 
pressure is indicated by the dashed line (jBertoldi McKed ll992), 
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Fig. 11. — Mass spectrum of the H CO + cores in the p Ophiuchi main cloud. The error bars 
correspond to yN counting statistics, where N is the number of cores in each mass bin. The solid 
lines are the least square fits in the interval of 1M© < M L te < 7M (AN/ AM oc M^ 3 ) and 
M LTE > 7M (AN/ AM oc M^). 
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Fig. 12. — Relationship between the surface term, S, and the gravitational term, W, in the virial 
equation. They are normalized to the internal kinetic energy term, U [Note that W/U = — (ovir) -1 ]- 
The solid line indicates the virial equilibrium, U + W + S = 0. The dashed line indicates the line 
at which W = S. For the cores that lie below the solid line, the value of U + W + S is negative 
and thus expected to be bound. All others are unbound and expected to disperse away, if they do 
not gain more mass through accretion and/or merging with other cores, or reduce internal support 
through turbulence dissipation. The filled squares and open circles are the same as those in Figure 

El 
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Fig. 13. — Line width-radius relation {upper panel) and mass-radius relation (lower panel) of the 
H 13 CO + cores toward the p Ophiuchi main cloud (open squares) and toward the Orion A molecular 
cloud (crosses). The p Oph cores identified from the smoothed data are indicated by filled circles, 
for comparison. The solid line in the upper panel indicates the best-fit power-law for the p Oph 
cores identified from both the original and smoothed data. The dotted and dashed lines indicate 
the Larson and the Heyer & Brunt relations, respectively. 
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Fig. 14. — Virial ratio-mass relation of the p Oph cores and the Orion A cores. For Orion A, 
the fractional abundances of H 13 CO + are set to 4.8 x 1CT 11 (adopted by Ikeda et al. 2007) and 
1.72 x 10~ n (adopted for p Oph in the present paper) in panels (a) and (b), respectively. The 
open squares, crosses, and filled circles are the same as those of Figure [13] T he dashed lines are 
the sa me as that of Fig. [TUl Note that the Jeans mass given by eq. [2.13] of iBertoldi fc McKee 
(| 19921 ) is almost the same for both the original and smoothed p Oph cores. 
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1 OT on o 

16 27 29.2 


-24 41 37.5 


3.67 


0.79 


86.1 


0.0522 


1.13 


0.525 


2.27 


5.37 


2.37 


0.66 


\r 
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52 


16 27 30.7 


-24 40 35.7 


4.32 


1.38 


84.0 


0.0509 


1.62 


0.316 


2.53 


3.37 


1.33 


0.79 
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53 


16 27 33.3 


-24 26 11.6 


3.93 


1.63 


114.2 


0.0692 


1.46 


0.563 


8.99 


7.72 


0.86 


1.12 
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16 27 33.3 
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4.58 


1.97 


95.6 


0.0579 


1.01 


0.361 


4.96 


4.21 


0.85 


1.06 
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3.80 


0.76 


92.1 


0.0558 


1.12 


0.675 


2.01 


7.85 


3.90 


0.48 
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16 27 36.5 
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3.41 


0.51 


76.6 


0.0465 


2.33 


0.603 


0.90 


5.64 


6.29 


0.37 


N 


57 


16 27 38.3 


-24 41 16.4 


3.54 


0.53 


57.7 


0.0350 


1.61 


0.331 


0.47 


2.39 


5.10 


0.45 


N 


58 


16 27 42.9 


-24 42 38.3 


3.80 


1.06 


76.8 


0.0466 


1.51 


0.444 


3.19 


4.03 


1.26 


1.31 


Y 
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16 27 46.0 


-24 44 41.6 


3.67 


1.31 


80.8 


0.0490 


1.53 


0.302 


2.97 


3.16 


1.06 


1.05 


Y 


60 


16 27 51.9 


-24 42 58.3 


3.93 


1.02 


75.1 


0.0455 


1.02 


0.341 


1.70 


3.17 


1.87 


0.75 


N 



Table 1 — Continued 



1 1 J 


R A 




VhSB, 


-^A,peak 


D 

J (-core 


p 

1 ^corc 


Aspect 


G^corc 


JWLTE 


M ■ 
1VI V1I 


tVvir 


71 






(J2000.0) 


(J2000.0) 


(km s^ 1 ) 


(K) 


(arcsec) 


(PC) 


Ratio 


(km s" 1 ) 


(M Q ) 


(M e ) 




(10 5 cm- 3 ) 




61 


16 27 57.7 


-24 35 4.8 


4.19 


0.97 


56.9 


0.0345 


1.44 


0.398 


1.31 


2.71 


2.06 


1.32 


N 


62 


16 27 59.2 


-24 34 23.6 


4.19 


0.75 


41.0 


0.0248 


1.01 


0.353 


0.64 


1.78 


2.80 


1.72 


N 


63 


16 28 0.6 


-24 33 1.2 


4.32 


1.34 


73.0 


0.0442 


1.05 


0.402 


2.43 


3.51 


1.44 


1.16 


N 


64 


16 28 3.7 


-24 34 43.8 


4.19 


0.75 


48.5 


0.0294 


1.68 


0.238 


0.60 


1.68 


2.79 


0.99 


N 


65 


16 28 6.7 


-24 34 23.1 


4.32 


0.71 


62.9 


0.0381 


1.14 


0.350 


0.91 


2.71 


2.97 


0.68 


N 


66 


16 28 15.8 


-24 35 44.7 


4.19 


1.06 


69.5 


0.0421 


1.04 


0.306 


1.38 


2.74 


1.98 


0.77 


N 


67 


16 28 20.4 


-24 36 25.6 


4.19 


1.60 


52.6 


0.0319 


1.11 


0.324 


1.94 


2.15 


1.11 


2.48 


Y 


68 


16 28 26.4 


-24 36 4.6 


4.19 


1.20 


97.4 


0.0590 


1.34 


0.271 


3.42 


3.58 


1.05 


0.69 


N 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. Cores 
No. 1 through 17 and 19 are located in the Oph A region and the temperature is set to be 18K. For the other cores, T = 12K. The typical 
uncertainties of 7? CO re and dv CO rc are, respectively, about 0.01 pc, derived from the uncertainty in the estimation of the core projected area, and 
the velocity resolution of 0.13 km s _1 . 

a Y and N mean, respectively, the cores associated with and without YSOs that are classified as either Class I or Flat Spectrum objects based 
on the Spitzer observations. 
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Table 2. Summary of the Physical Properties of the H 13 CO + cores 



Property 


Minimum 


Maximum 


Mean a 


Median 


^corc (pc) 


0.0221 


0.0692 


0.0450 ± 0.0113 


0.0451 


dv COTC (km s -1 ) 


0.186 


0.765 


0.488 ± 0.142 


0.474 


M LTE (M ) 


0.39 


21.79 


3.35 ± 3.57 


2.42 


M vir /M LTE 


0.368 


7.76 


2.37 ± 1.57 


1.89 


n (cm -3 ) 


3.70 x 10 4 


5.77 x 10 5 


(1.35 ±0.96) x 10 5 


1.35 x 10 5 


Aspect Ratio 


1.01 


2.33 


1.34± 0.33 


1.22 



With standard deviation 



Table 3. Summary of the Physical Properties of the 850/um cores 



Property 


Minimum 


Maximum 


Mean a 


Median 


Rcore (PC) 


0.00494 


0.0601 


0.0204 ± 0.0118 


0.0174 


Mlte (M ) 


0.0868 


8.25 


1.42 ± 2.13 


0.553 


n (cm -3 ) 


1.41 x 10 5 


3.20 x 10 6 


(6.53 ± 6.29) x 10 5 


4.37 x 10 5 



1 With standard deviation 



Not e. — The uniform temperature of T = 15 K is assumed for all the cores in Ijprgensen et al 



(2008). 
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Table 4. Summary of the Physical Properties of the H 13 CO + cores with and without YSOs 



Property 


Minimum 


Maximum 


Mean a 


Median 


Cores without YSOs 


#corc (PC) 


0.0221 


0.0612 


0.0403 ± 0.0105 


0.0398 


dv corc (km s _1 ) 


0.186 


0.722 


0.431 ± 0.132 


0.421 


M LTE (M ) 


0.469 


6.89 


1.54 ± 1.69 


1.54 


M vir /M LTE 


0.74 


7.76 


2.63 ± 1.65 


2.05 


n (cm -3 ) 


3.70 x 10 4 


3.31 x 10 5 


(1.29 ±0.66) x 10 5 


1.30 x 10 5 


Aspect Ratio 


1.01 


2.33 


1.45± 0.417 


1.36 


Cores with YSOs 


-Rcorc (pc) 


0.0295 


0.0692 


0.0497 ± 0.0102 


0.0492 


C orc (km s _1 ) 


0.302 


0.765 


0.544 ± 0.130 


0.545 


M LTE (M ) 


0.394 


21.8 


4.55 ± 4.48 


3.55 


M vir /M LTE 


0.37 


5.54 


2.11 ± 1.46 


1.53 


n (cm -3 ) 


0.449 x 10 5 


5.77 x 10 5 


(1.40 ± 1.20) x 10 5 


1.12 x 10 5 


Aspect Ratio 


1.02 


1.62 


1.22± 0.155 


1.19 



a With standard deviation 



